Gutzwiller-projected wave functions for the pseudogap state of underdoped high-temperature 

superconductors 
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Recent experiments strongly suggest tiiat a Fermi surface reconstruction and multiple Fermi pockets are 
important common features of the underdoped high-temperature cuprate superconductors. A related theoretical 
work [Phys. Rev. B 79, 134512 (2009)] has demonstrated that a number of hallmark phenomena observed in 
the underdoped cuprates appear naturally in the scenario of a paired electron pocket co-existing with unpaired 
hole pockets. We propose Gutzwiller-projected wave-functions to describe this two-fluid state as well as two 
competing states in its vicinity. It is argued that a pseudogap state constructed from these wave-functions may 
be selected by energetics at finite temperatures due to spin fluctuations. 



Shortly after the discovery of the first high-temperature 
cuprate superconductor, Anderson recognized' strong elec- 
tronic correlations as the key to understanding a mechanism 
for superconductivity in these compounds and proposed a 
Gutzwiller projected BCS wave-function as an approximate 
ground state. These RVB ideas have been pursued further 
in a number of follow-up works,!^ which included analyt- 
ical methods of renormalized mean-field theory,"* numerical 
Monte-Carlo variational analyses,'3Si as well as effective field 
theory approaches where non-linear RVB constraints were im- 
plemented via an auxiliary gauge field.^ Another key aspect of 
the cuprate phase diagram, the presence of an antiferromag- 
netic (AF) Neel order at low dopings, was emphasized later 
in the framework of the spin-fermion (SF) model,** '^ where 
strong AF fluctuations in the neighboring metallic normal 
state were proposed as the pairing glue. 

The recent experimental discovery'*" of small Fermi 
pockets in the underdoped cuprates has initiated a renewed 
interest in t he p hysics of these amazing materials and led to 
new ideas. '^'^ In particular, the original breakthrough paper'" 
suggested that the Fermi pockets observed via quantum oscil- 
lations might be electron rather than hole pockets. A later the- 
oretical worlP^ showed that the "nodal-antinodal dichotomy" 
observed in experiment arises naturally within the picture of a 
strongly-paired (but uncondensed) electron pocket in the anti- 
nodal region and unpaired hole-pockets in the nodal regions. 
The approach of Ref. Il20l is to drive a magnetically ordered 
AF metal state (with electron and hole excitations naturally 
present due to a Brillouin zone folding) into a critical region 
via strong AF fluctuations, which were shown to give way to 
a strong i-wave pairing in the electron pocket and a weaker 
p-wave pairing in the hole-pockets. The conventional d-wave 
symmetry is restored via a Brillouin zone unfolding. 

Although spin fluctuations are the basic pairing glue in the 
SF model on both the normal Fermi-liquid side^ and the AF 
metal side,^*^^' they do not induce pair-breaking on the anti- 
nodal electrons in the latter model.^" The SF model in an AF 
metal is effectively of Ising type^^ and the s-wave electron 
pairing is the strongest right in the critical point in a sharp 
contrast to the normal Fermi-liquid result, where the pair- 
breaking dominates at criticality. This suggests that it may be 
energetically favorable to retain the AF coiTelations to a mini- 
mum necessary to give rise to a folding and induce the pairing 
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FIG. 1: (a): Separated nodal hole pockets (solid red line) and antin- 
odal electron pockets (dashed blue line) in a doped antiferromagnet. 
(b): In the paramagnetic phase the hole and electron pockets touch 
each other. The thin dashed line marks the magnetic Brillouin zone 



at criticality. The energy gain due to the BCS pairing can com- 
pensate the presence of otherwise unwanted AF correlations, 
such that the charge sector stabilizes a finite-temperature spin- 
liquid state. 

In this paper, we attempt to bring together the SF model 
and the RVB theory by proposing a set of Gutzwiller projected 
wave-functions to construct the pseudogap state of Ref. 1201 . 
We note that our approach below is distinct from other similar 
variational wave-function methods used to study competition 
of antiferromagnetism and superconductivity,^-' ^^ as we treat 
the nodal and antinodal quasiparticles differently. Our starting 
point is the f - f ' - /-model 
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where c\g. creates an electron of spin cr on site r and Sr - 
5 2 Cro-To-o-'Cro-' is the corresponding spin operator, with Ta-a' 

being the Pauli matrices. Here, t and t' are respectively the 
nearest and next-to-nearest neighbor hopping amplitudes, and 
J - 4t^/U is the nearest-neighbor AF superexchange cou- 
pling, which is related to the on-site Hubbard repulsion, U. 
The t - t' - J model has been widely used to study the effects 
of strong correlations in the superconducting cuprates. On the 
other hand, this is also the minimal model which can describe 
antiferromagnetic coiTelations due to the superexchange term. 
A mean-field Hamiltonian that gives rise to an AF metal 



state can be constructed from the t-t' - 7-model ([TJ as follows 

'HaFM = - 2_i fr,r'cio.Cr'o- - -/z / ^ 
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where z = 4 is the co-ordination number, (-1)^*^"^' = ±1 if r e 
Saib, with S AjB denoting the A/B-sublattices. In Eq. (pi, m is 
the modulus and the unit vector n is the direction of the Neel 
order parameter, which we choose to be in the z direction. 
Using the Fourier transforms fiko- - ^res^ Cro-e "'' and 

by a- - YjteSs ^ro-s "''i ^nd the four-component operator *?J^ - 
(ak,t>^k,T>^k,i,^k,i)' the Hamiltonian Ob reads 

"Hafm = Yj *^I[«*)'''- X 4 + £1 (k)/, X &, + 0f, X (r,]>?u, (3) 



where e(k) = -2n cos A:^ -H cos A:^ j and ei(k) = 
At' cos kx cos A;j, are the inter-sublattice and intra-sublattice 
dispersions, & and f are the Pauli matrices acting respectively 
on the sublattice and spin indices, and I^ and /o- are the 
2x2 identity matrices. Hamiltonian (|3]l is diagonalized 



by a Bogoliubov rotation and yields two eigen-modes, 
^iCk) = fif/ftCk) = ei(k) ± Ve^(k) + (j? with = 7zm/2, 
where "e" and "h" correspond to electron- and hole-branches 
located primarily in the antinodal and nodal regions of the 
magnetic Brillouin zone (MBZ) respectively, as shown in 
Fig.fTI The electron and hole operators are given by 



■k+ 



= [a^.Ia + '^ko-y] ,. and ,. = [p^I^ + '^kO-^] ,. 

"k+'' ^^kT*^ ^ k-'' ^\i, 

(4) 
where a^ = ^ [i _ 0/ ^sH^i) + f] and \a\,\^ + ^kl^ = 1- Note 
that the shapes of the electron and hole pockets depend on the 
details of the lattice hoppings via the ratio tjt' , and we assume 
that it is such that both pockets appear at low dopings.E^I^^I 

We now use operators (J4|i to express variational wave- 
functions describing a paired electron pocket and competing 
metallic and superconducting states. We denote these states 

by m, Ae,Ah), where m is the magnetization, and Ag/h is the 
electron/hole variational gap parameter The simplest such 
wave-function describes an unpaired AF metal 
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where the Gutzwiller projector, Pq - HrC^ ~ «rT«ri), enforces the non-linear no-double-occupancy constraint, ^(./^(k) are 
the spectra given above Eq. (Hli, and yu is the chemical potential. Here and below the overline implies that the corresponding 
parameter {e.g., magnetization, m, in Eq.|5]l is calculated self-consistently by minimizing the energy. 

Another possible wave-function, argued to be of relevance to the pseudogap state, is a Gutzwiller-projected product of wave- 
functions for the paired electron pocket and unpaired hole pockets: 



|pPEP) = |m,Ae,0)=!PG W [U\ 
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where f/^ and V^'' are the usual Bogoliubov amplitudes for a pairing in the electron pocket. I.e., \U^\ + V^"^ - 1 and 
|-y{c)j2 _ 1 N _ [£'^(k) - ji\l ^[EeQ/i) - lA^ + A^(k)|. The variational parameter, Ag, corresponds to the i-wave symmetry of the 
electron pairing gap in the folded Brillouin zone. But it is an anisotropic i-wave in the MBZ, if we identify the antinodal regions 
connected by a reciprocal lattice vector Q = {n,n). Note that a non-zero Ag before projection does not guarantee a long-range 
order, which is strongly suppressed by the non-linear Gutzwiller projection.13 

For completeness, we also write down the variational wave-function for an exotic fully-paired state 
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where f/j^ and V^ are the Bogoliubov parameters in the hole 
pockets. As noted in Ref. [20], the hole-pairing must have a p- 
wave character in order to reproduce the known (i-wave sym- 
metry in the unfolded pictm^e. However,only in the absence of 
a long-range Neel AF order and if both the electron and hole 
gaps are non-zero, the state reduces to the more conventional 
projected li-wave BCS state expressed in terms of the physical 



electrons | PBCS) = Pg Ukesz [Uk + Vkt'l/l^^^] \ VAC). 

We now focus on states (|5]l and (|6| to calculate their 
energy and determine the optimal variational parameters in 
the fram ework of Gutzwiller renormalized mean-field theory 
(RMFT).EEElThe basic assumption of RMFT is that the ex- 
pectation value of any operator, O in the projected state is pro- 
portional to that in the corresponding unprojected state, i.e.. 



(PiA 1^1 PiA) = {^ |!Pg^!Pg| iA) ~ 0o('/' |<5| (A)- The pro- 
portionality factor, gr^), is called the Gutz wilier factor. The first 
step in building RMFT is to determine the Gutzwiller factors 
for all operators in the tJ Hamiltonian (fill. 

We first consider the relation between the on-site den- 
sities, Hr.o-, in the projected states and that in the unpro- 
jected states, n|.^. In an AF-ordered state, we have n^'^ - 

ng_^ - («o/2)(l + crm), where m is the staggered magne- 
tization in the unprojected states. Note that the unprojected 
AF metal state is an eigenstate of the total particle number, 
hence ( PAFM I A^ I FAFM ) = ( AFM | A^| AFM). There- 
fore, we have in this case 2a-«r,o- - Ho-'^r.cr, which leads 



to Hr.o- = ««r,(T 1 ~ "r-cT / " ~ 2«i._o.«j. _^ , where the "bare 



density," n, is related to the doping level, x, simply as n = 1 -x. 
Consider now the hopping terms cl^Cr^a- In a projected 
state, this term contributes in a configuration which has hole 
at site r and a spin cr at site r' and the action of the operator 
is to reverse this configuration r <-> r'. The probabilities for 
such a process to occur in a projected//unprojected states are 
oc [rir^^ (1 - nJInr^a) «r',tr (1 - "r'/Zwr',^)]'^^- Their ratio de- 
termines the Gutzwiller factor for hoppings: 

9^" -ii-n) f..«rw[e(i-ie)<v(i-'C)]- 

Using the relation between nr,o- and rirl derived above, 
we obtain the following Gutzwiller factors for inter-sublattice 



hopping gf^ = gi - Ix/rf^ (x) (which is spin-independent), 

and the spin-dependent inra-sublattice hopping g, - 

gf'^ = giO) = ^4^'i^)l [rit\M[\^)\ and g^^'^ = gf'^ ^ 

g-i — Ixif^ {x)l\ri.^ \x)ti^ \x)\, where we introduced the 

functions rj'^ \x) - \ + x ±m}^(\ - x) for brevity. The spin- 
flip renormalization factor can be obtained similarly and reads 

gfi = 4 J]2 \x) ■ Interestingly in the limit of a nearly-perfect 
Neel spin polarization, m — > 1, at low dopings (x <K 1), we 
find gi ~ X, g, ~ g2 ~ I, and g^ ~ x^. This is consistent with 
the idea by Kane et al?^ that holes in an AF background are 
coherent due to intra-sublattice hopping. We note that in eval- 
uating the Gutzwiller factors, we have only taken into account 
on-site correlations due to projection and neglected non-local 
correlations of the wavefunctions.^^ 

Using these Gutzwiller factors, we calculate the energy of 
a generic variational projected state as follows: 
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FIG. 2: (Color online) Neel order parameter in the Gutzwiller- 
projected antiferromagnetic metal ipl as a function of doping for 
three different sets of hopping parameters. In general, a larger t' 
sustains magnetization over a wider range of doping levels, as the 
energy gained from intra-sublattice hopping favors the AF state. 

P(k,k') = [370,/(8O]e(k - k'), and 2(k,k') = 



81, the 



keMBZ,(T;v=± 
k.k'eMBZ.tr 
k.k'eMBZ 
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where 5+(k) = ei(k)(g2 + gs) /2, 5_(k) 
2<7ie(k)akj6k + ei (k) (/^^ - ck2) (^2 - ^3) /2, /(k,k') = 
[JgJ{2t)] [(3/2)e(k - k')ak;Skak'y8k' -t{pl- ff^) [pI - 4)], 



-[7g,/(80]£(k - k')(yS2 - a'^ifil - 4). In Eq 
functions n(,(k, cr) and n/,(k, cr) are the occupation numbers 
of the electrons and holes in the corresponding unprojected 
states. E.g., in the AF metal state, they are simply the 
Fermi-Dirac distributions. 

We first focus on the unpaired AF metal state J5]), where 
VvSv. — 0. Minimizing the energy functional ([8k we ob- 
tain the staggered magnetization as a function of doping. The 
magnetization goes down with doping with the precise slope 
and the critical doping at which the magnetization vanishes 
being dependent on the value of t' jt, see Fig.l2] 

We next look at the paired-electron-pocket wave- 
function (|6]l. An energetics analysis shows the electron 
gap, Ae, increasing monotonically with underdoping (c./, 
Ref fSl) and indicates that small but finite values of m lead 
to higher energies, hence ih — > -(-0 is the optimal parameter 
at the mean-field level. Fig. 3] compares the 7 = energies 
of the projected AF metal (|5 1 and the paired-electron-pocket 
state (|6]l for a specific reasonable choice of the hopping 
parameters. Interestingly, the two energies cross, indicating 
that in the corresponding window of doping levels, a par- 
tially paired state is energetically preferable. However, the 
precise value and even the existence of this window depends 
strongly on the details of quasiparticle dispersion and hence 
is non-universal and is strongly susceptible to fluctuation 
renormalizations . 

Finally, we study the wave-function ^\ with both electron 
and hole pockets paired in the absence of AF order. The state 
reduces to the projected li-wave paired state written in terms 
of the physical electron operators and with (A*^ H- A'') being an 
effective li-wave gap. As expected, this fully paired state has 
the lowest energy at zero temperature and yields the projected 
£/-wave superconductor as the ground state !^ 

However, this does not negate the consideration of "pseu- 
dogap" wave-function (J6]l above, since the latter intends to de- 
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FIG. 3: (Color online) Displayed are the energy of the paired- 
electron-pocket state (ISl and the energy of the projected antiferro- 
magnetic metal Jsl with optimal magnetization and for the hopping 
parameters: t = 2J and t' = -t/4. Note that at low but finite dop- 
ing, the two energy curves cross indicating that the state J6l becomes 
energetically preferable to the fully unpaired AF metal pT 



scribe finite-temperature phenomena, which the pseudogap is. 
Let us invoke here the arguments of Ref. 1201 . which showed 
that spin fluctuations have a detrimental effect on the p-wave 
pairing of the holes. Note that the conventional projected 
wave-function approach greatly reduces the Hilbert space to 
that of trial functions, which do not include excited states. 
Hence, it can not properly capture thermal fluctuations and 
spatial inhomogeneities. However, one can construct a "vari- 
ational density matrix" within this trial Hilbert subspace 

d^md'-Ase -r m, Ae,0)(m, Ae,0 (9) 

to model the impact of (uniform) magnetization and pairing 



fluctuations. Here, T is temperature and Tim, Ae) - amin^ -i- 
/3„,m'* + Ue \K?' + Pe \Aet + J \K\^ Hi^ is a Landau-type func- 
tional with a's, jS's, and y being the variational parameters de- 
termined by minimizing the energy, £[;F] = Tr rKtjp /Trp, 
where the trace is over the truncated Hilbert space. Note that 
in general, a non-zero fluctuating Ah is also allowed in (J9]l. 
However, such a state does not reduce to the projected t/-wave 
superconductor, and therefore the fluctuating gaps Ag and Ah 
are described by two independent sets of variational param- 
eters. If AF fluctuations are strong enough, the pseudogap 
state (J9|l (with Ah being strongly suppressed) is selected by 
energetics. We note here that a proper quantitative analysis 
of the energetics should include spatial dispersion for fluctua- 
tions, restoring a larger Hilbert space and bringing us back to 
a field theory-ISSl 

In conclusion, we have proposed a set of Gutzwiller pro- 
jected wave-functions (|6]) and a related "variational denisty 
matrix" (J9]l, where unpaired nodal holes co-exist with paired 
antinodal electrons, to describe the pseudogap state of the un- 
derdoped high-temperature cuprate superconductors. Finally, 
we mention that a smoking-gun experiment that could confirm 
the pseudo-gap state proposed in Ref. Il20ll and here would vi- 
sualize the electron pocket in the underdoped region by sup- 
pressing the strong pairing there. Apart from the transport 
measurements in ultra-high magnetic fields,'" '^ this could 
also be achieved by performing ARPES measurements in the 
non-linear regime of high currents (exceeding the critical cur- 
rent for the electron pocket) that should destroy a strong s- 
wave pairing near the antinodes and open up the hidden elec- 
tron pocket. 

Acknowledgements: VG is grateful to Professor Anderson 
for an illuminating discussion. The authors are indebted to 
Professor Sachdev for reading the manuscript and providing 
valuable comments. This research was supported by the De- 
partment of Energy (VG) and DARPA-QuEST (RS). 



P. W. Anderson, Science 235, 1196 (1987). 

G. Kotfiar and J. Liu, Phys. Rev. B 38, 5142 (1988). 

B. Edegger, V. N. Muthukumar, and C. Gros, Adv. Phys. 56, 927 

(2007). 

F C. Zhang, et al, Supercond. Sci. Tech 1, 36 (1988). 

A. Paramekanti, M. Randeria, and N. Trivedi, Phys. Rev. B 70, 

054504 (2004). 

P W. Anderson, etal.J. Phys.: Condens. Matter 16, R755 (2004). 

P A. Lee, N. Nagaosa, and X. G. Wen, Rev. Mod. Phys. 78, 17 

(2006). 

V. Barzykin and D. Pines, Phys. Rev. B 52, 13585 (1995). 

A. Abanov, A. V. Chubukov, and J. Schmalian, Adv. Phys. 52, 1 19 

(2003). 

N. Doiron-Leyraud, et al. Nature 447, 565 (2007). 

D. LeBoeuf, et ai. Nature 450, 533 (2007). 

S. E. Sebastian, et ai. Nature 454, 200 (2008). 

A. J. Minis and M. R. Norman, Phys. Rev. B 76, 220503(R) 

(2007). 
'" K. Sun, etal., Phys. Rev. B 78, 085124 (2008). 
15 C. Varma, Phys. Rev. B 79, 085110 (2009). 



T. Senthil and P A. Lee, Phys. Rev. B 79, 2451 16 (2009). 

S. Sachdev, Physica Status Solidi B 247, 537 (2010). 

X. Jia, P Goswami, and S. Chakravarty, Phys. Rev. B 80, 134503 

(2009). 

T. Das, R. S. Markiewicz, and A. Bansil, Phys. Rev. B 81, 184515 

(2010). 

V. Galitski and S. Sachdev, Phys. Rev. B 79, 134512 (2009). 

E. G. Moon and S. Sachdev, Phys. Rev B 80, 035117 (2009). 

M. Dzero and L. P Gor'kov, Phys. Rev. B 69, 092501 (2004). 

T. Giamarchi and C. Lhuillier, Phys. Rev. B 43, 12943 (1991). 

S. Pathak, et al, Phys. Rev. Lett. 102, 027002 (2009). 

M. Ogata and H. Fukuyama, Rep. Prog. Phys. 71, 036501 (2008). 

B. L. Altshuler etal., Europhysics Lett. 41, 401 (1998). 
A. V. Chubukov and D. Morr, Phys. Rep. 288, 355 (1997). 

C. L. Kane, P A. Lee, and N. Read, Phys. Rev. B 39, 6880 (1989). 
M. Ogata and A. Himeda, J. Phys. Soc. Jpn. 72, 374 (2003). 

V. B. Geshkenbein, L. B. loffe, and A. I. Larkin, Phys. Rev. B 55, 
3173(1997). 



